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ABSTRACT 


This report presents shadow data and heat input data for the S-49 
EGO. It gives times for which each experiment is in the shadow of 
either the Earth, the satellite's main box or the solar array. It also 
gives the heat inputs to the experiments as a function of time from 
launch for one complete orbit. The heat inputs include direct solar 
radiation, reflected solar radiation and Earth emitted radiation. 

The shadow data are presented for a launch date of May 15, 1964 
at 13.0 hours U.T. The heat input data are for a launch time of May 
15, 1964 at 8.0 hours U.T. The injection elements are from Reference 
2. The satellite's geometry is from References 3 and 4. 
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I. INTRODUCTION 


The Orbiting Geophysical Observatory (OGO) consists of a main 
box, the solar array and the orbital plane experiment package (OPEP). 
Experiments may be appended to any of the three main component parts 
either directly or on booms. As many as fifty experiment packages 
may be carried on the OGO spacecraft. (See Figure 1). The solar 
array and the OPEP may rotate relative to the main box. The main box 
and the solar array are considered as casting shadows and carrying ex- 
periments, the OPEP is considered only as carrying experiments. A 
method for calculating the radiant heat inputs to each experiment is 
also presented in this report. These inputs include direct solar, re- 
flected solar and Earth emitted heat flux. If an experiment is shielded 
from the Sun by the Earth, the main box or the solar array, the solar 
input to that particular experiment is set to zero. If an experiment is 
in the eclipse of the Earth the reflected solar (as well as the direct 
solar) heat input is set equal to zero. However, the effects of shadow- 
ing of an experiment from the Earth by the main box or the solar array 
is not considered in the computation of the solar reflected or Earth 
emitted heat inputs to that particular experiment. Shadow data are 
presented for the S-49 4^GO) for a launch date of May 15, 1964 at 13.0 
hours U.T. The heat input data are for a launch time of May 15, 1964 
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at 8.0 hours U.T. The injection elements are from Reference 2. The 
S-49 geometry is from References 3 and 4. 
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II. NOTATION 


A Cross sectional area. 

a The fraction pf the solar constant a which is re- 
flected away from the Earth (albedo). 

a E The equatorial radius of the Earth. 

f « (7i> 7 2 > 73 . r ) The integrand of the integral for q na (t). In 

equation (47) it is defined on page 38. 

f E ( y 2 , y 3 , r) The integrand of the integral for q nE (t). In 
equation (60) it is defined on page 44. 

f(i) Incidence factor defined by equation (45). 

H = r VE x r yE ; the angular momentum vector per unit 
mass. 

I(t) Shadow indicator at time t ; 0 in shadow, 1 not in 
shadow. 

i Angle of incidence, 
i, J, k Orthogonal unit vectors. 
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= (a,b,c) ; the 3-vector representing the line 
connecting the projections of points j and j + 1 onto 
a test plane, as defined by equations (20, 21), and 
appendix C. 

s p." +l - p." ; as used in appendices C and D. 

Outward normal to an elemental area AA 7 of the 

i m 

Earth’s surface; it is defined in equation (49). 

Unit vector normal to an experiment surface. 

Augmented position vector (x/', y.", l) . 

A position vector ( x , y/ # ) defined by equation (17). 

Augmented position vector (x n " , y n " , 1 ) . 

A position vector (x n ", y n ") defined by equation (18). 

A position vector used in appendix B. 

The heat flux incident upon a surface of the n th ex- 
periment at time t , which is emitted by the Earth. 

It is given by equations (60) and (62). 
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q nE The accumulated heat flux q nE over a time interval; 
it is defined by equation (63). 

q n , (t) The heat flux incident upon a surface of the n th ex- 
periment at time t , which originates at the Sun and 
reflects from the Earth. It is defined by equations 
(47) and (58). 

q na The accumulated heat flux q na (t) over a time 
interval; it is defined by equation (59). 

q ns (t). The heat flux incident upon a surface of the n th ex- 
periment at time t, which radiates directly from 
the Sun. It is defined by equation (43). 

q ng The accumulated heat flux q ns over a time interval; 
it is defined by equation (46). 

R Augmented position vector (x, y , z, 1 ), (a 4-vector). 

r A position vector (x, y, z ). 

r yE The vector which points from the Earth to the 
vehicle. 

r yE The velocity of the vehicle relative to the Earth. 
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r lm The vector which points from the satellite to the 
Earth's elemental surface area AA im . It is given 
by equation (51). 

r sy The vector which points from the vehicle to the 
Sun. 

T sim The vector which points from the Earth's elemental 
surface area AA 7m to the Sun. 

fr m 

r gE The vector which points from the Earth to the Sun. 

S i s B i> C i ,D i ) ; a 4-vector representing a face 

defined in equation (3). 

s Defined by equation (10). 

t Time. 


t 


o 


u 


u 


• Y Imn> Z 


k l mn » 1 Imn 


l mn 


Reference time. 


= (1, 1, 1); a vector. 


Defined by equation (22). 


Vectors ( x„ x m , x„) , (y,, y m . y„ ) and (zj, * m , z„) , 
respectively. 
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x, y, z The coordinate system associated with i, j , k , 
respectively. 

a The angle corresponding to the longitude in the 

plane normal to r yE , the zero line being the projec- 
tion of i fi on this plane. (See Figure 5). 

a The value of a associated with the elemental area 

m 

AA i m ; defined by equation (50). 

ft The angle between r yE and the tangent to the Earth 
from the satellite; (equation 48). 

y x The angle of incidence of the Sun’s rays on an 
elemental area of the Earth's surface. 

y lim The value of y x associated with the elemental area 
AA i m . It is the angle between the vectors N^ m and 
r slm (Figure 5). It appears in equation (52). 

y 2 The angle between the radius vector which points 
from an elemental area of the Earth's surface to 
the vehicle and the normal to that area. 
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y 2 l m The value of y 2 associated with the elemental area 
AAj m . It is the angle between the vectors N [m and 
(- 7 lm ) (Figure 5). It appears in equation (54). 

y 3 The angle between the vector which points from 
the vehicle to an elemental area of the Earth's 
surface and the unit outward normal to the ex- 
periment surface. 

y 3lm The value of y 3 associated with AA Im . It is the 

angle between the vectors n and r (Figure 5). It 
appears in equation (55). 

AA Jm An incremental surface area of the Earth (Figure 5); 
it is defined in equation (57). 

d The angle between the vector r yE and the radius 
vector from the center of the Earth to the surface 
element in question. 

The value of 6 associated with the elemental area 
AA Jm . It is the angle between and r yE . Equa- 
tion (50). 


8 


A The emittance per unit area of the surface of the 
Earth. 

c r The solar constant 

<f> The solar array angle; defined by equation (35). 
(Figure 1). 

V The OPEP angle; defined by equation (42). (Figure 1). 



Superscripts and Subscripts 

Denotes rotated coordinates in equations (4), (5), 
(7), (8) and (9). Also denotes a change of origin 
location in equations (14) and (15). 

Denotes vectors in the test plane. 

Box. 

Referring to the n th experiment. 

Solar array. 

Referring to a point on the solar array axis. 
OPEP. 


Referring to a point on the OPEP axis. 



Vector Operations 


Ax B The cross product of the sectors A andB . 


A • B The scalar product of vectors A andB . 


^P ^P 

W = T ^ +T + k faT ; the gradient of the scalar 


function F . 
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HI. ASSUMPTIONS 


A. Initial Conditions 

The initial conditions are assumed to be the same as the Agena 
burnout conditions (Reference 2) namely: 

| 

geocentric latitude = -20.28 degrees 

terrestrial longitude = 1 1 1.902 degrees 

geocentric height = 281621.6 meters 

the speed of the vehicle = 10645.27 meters per second 

azimuth = 66.355 degrees 

flight path angle = 1.47148 degrees 

B. Restraints on Launch Time 

The results presented in this report are for launch times for which 
perigee does not go below its initial value for one year. See Reference 
5 for the results of a detailed launch window study. 

C. S-49 Spacecraft Geometry 

The assumed geometry, presented in this section, is given in the 
box coordinate system. The center of gravity (c.g.) of the spacecraft 
was arbitrarily chosen as the origin. The coordinates of the solar 
paddles and the OPEP are given in a standard position relative to the 
box (solar array angle $ = 0 and OPEP angle m = 0) together with the 
location of their axes of rotation. 
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1 . Shadowing Surfaces 


a. Main Box — The box coordinates of the vertices (corners) 
of the main box are a follows: (See Figure 2). 


Point 

X 

inches 

y 

inches 

z 

inches 

1 

-14.607 

23.5 

14.983 

2 

-14.607 

-43.5 

14.983 

3 

14.607 

-43.5 

14.983 

4 

14.607 

23.5 

14.983 

5 

-14.607 

23.5 

-14.983 

6 

-14.607 

-43.5 

-14.983 

7 

14.607 

-43.5 

-14.983 

8 

14.607 

23.5 

-14.983 


b. Solar Array — The coordinates of the solar array are as 
follows: (See Figure 2). 


Point 

X 

inche s 

y 

inches 

z 

inches 

PI 

-117.841 

0.0 

-36.877 

P2 

-117.841 

0.0 

36.877 

P3 

- 39.541 

0.0 

36.877 

P4 

- 39.541 

0.0 

18.677 

P5 

- 26.491 

0.0 

18.677 

P6 

- 26.491 

0.0 

-18.677 


i 
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Point 

inches 

inches 

inches 

P7 

- 39.541 

0.0 

-18.677 

P8 

- 39.541 

0.0 

-36.877 

P9 

117.841 

0.0 

-36.877 

P10 

117.841 

0.0 

36.877 

Pll 

39.541 

0.0 

36.877 

P12 

39.541 

0.0 

18.677 

PI 3 

26.491 

0.0 

18.677 

P14 

26.491 

0.0 

-18.677 

PI 5 

39.541 

0.0 

-18.677 

PI 6 

39.541 

0.0 

-36.877 


The axis of rotation of the solar array passes through the point 

(V V z p„) 1 (°' °'°)- 

2. Experiments 

The coordinates of the experiments and their normals are as 
follows (See Figure 2). 

COORDINATES 


Experiment 

X 

y 

z 

Inches 

Inches 

Inchei 

E.P. 1 

- 74.2 

91.1 

12.3 

E.P. 2 

49.9 

100.1 

15.3 

E.P. 3 

- 12.7 

103.1 

- 4.2 
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Experiment 

X 

Inches 

y 

Inches 

z 

Inche s 

E.P. 4 

20.6 

- 94.4 

-19.2 

E.P. 5 

11.4 

288.1 

-25.2 

E.P. 6 

- 6.7 

-288.9 

-16.7 

SOEP 1 

-102.7 

3.4 

- .30 

SOEP 2 

102.4 

3.4 

- .30 

OPEP 1 

.2 

43.1 

45.8 

OPEP 2 

.2 

43.1 

-45.8 

axis of rotation of the 

OPEPS passes through the point 


, y T , z T \ - ( 0 , 43. 5» 

O O / 

o). 



Normals 

X 

y 

z 

E.P. 1 

-1.0 

0.0 

0.0 

E.P. 2 

0.0 

1.0 

0.0 

E.P. 3 

0.0 

1.0 

0.0 

E.P. 4 

0.0 

-1.0 

0.0 

E.P. 5 

0.0 

1.0 

0.0 

E.P. 6 

0.0 

-1.0 

0.0 

SOEP 1 

0.0 

1.0 

0.0 

SOEP 2 

0.0 

1.0 

0.0 

OPEP 1 

1.0 

0.0 

0.0 

OPEP 2 

1.0 

0.0 

0.0 


IV. MATHEMATICAL MODEL FOR SHADOW TESTING 


A. General Description of the Shadow Computations 

In order to carry out the shadow computations, the geometry of 
the satellite (the main box, the solar array, and the experiments) has 
to be described analytically. The configuration of the satellite is 
described in a coordinate system which is fixed in the main box. The 
shadowing of experiments by the main box and the solar array is de- 
termined by testing the shadow cast by each of their component faces. 
This may result in some duplicate testing, since part of the faces may 
overlap, but no special assumption need be made as to its shape. Any 
three dimensional object may be approximately described by a number 
of plane faces. 

The geometry of the satellite is input to the program by loading 
the coordinates of all its vertices and identifying the vertices belong- 
ing to each component face. Care must be taken that the polygon 
circumscribing each component face be convex. The occurrence of 
re-entering edges can and must be avoided by subdividing the cor- 
responding face into several convex polygons. The coordinates of 
points on the solar paddles and the OPEP are given in a standard 
position relative to the box (solar array angle $ = 0 and OPEP angle 
V = 0) together with the location of their axes. The program will then 
compute the rotation angle of the solar array and OPEP from their 
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attitude control logic and obtain the position of shadowing surfaces and 
experiments attached to paddles and OPEP. 

The program operates by transforming the coordinates of the Sun 
and the Earth into the box coordinate system by means of a transforma- 
tion matrix described in section IV-C. 

Fewer points are involved in this transformation than in the re- 
verse process. In this box coordinate system, the shadow cylinders of 
all the faces are then computed and the location of the experiments 
within or outside these cylinders is tested. If any experiment is found 
to be inside the shadow cylinder of any one surface, this experiment is 
in the shadow and no further testing is necessary for this experiment. 
Similarly, if the satellite is in the shadow of the earth, no shadow 
testing will be performed for any of the faces. 

B. An Experiment Shadowed From the Sun by a Polygonal Plane Face. 

1. The Equation of a Plane Formed by a Surface of the Main Box. 

The vertices of the box are denoted by r . ( i = 1 .2, • • • N). Out 
of these N points various sets of k points (k > 3) are given as vertices 
of the corresponding faces. The first three points of each such set are 
used to determine the equation of the plane in which the face lies. 

Care must be taken in the selection of these points, so that an 
accurate description of the surface is obtained from these points. In 
particular they must be chosen to lie on edges which are not too short 
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and not near parallel. This can always be achieved by arranging the 


points in proper order. If the first three points on the surface are 
denoted by r t , r m and r n the equation of the surface may be written in 
the form (see Appendix A). 


X 

y 

z 

l 

X Z 

Vl 

z l 

l 

X m 

y* 

Z m 

l 

X 

n 

y n 

Z 

n 

l 


( 1 ) 


or alternately by: 


Ai x + Bj y + Cj z + Dj - 0 (2) 

The four vector (Aj, B it C i , Dj is denoted by Sj and is taken as the 
analytic representation of the plane containing the i th face. The 
vectors r, r ; , F m and 7 n are in a coordinate system which is fixed to 
the main box. If the three -vector (xj, x m , x n ) is denoted by X lmn with 
the corresponding notation for Y Imn and Z lmn and U is defined as the 
vector U = (1, 1, 1) (see Appendix A); consequently. 

*i = 0-(v .„**>_) 

B, = 0 • (!,„ * X.J 

(3) 

c, = o • (x,.*Y,„) 

^Imn f^lmn * ^ Imnj 
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2. The Equation of a Plane Formed by a Surface of the Solar Array. 


r. and r. # represent position vectors of any solar array point in 
its zero and rotated position respectively. r po represents the coordinates 
of any point on the axis of rotation. The axis of rotation is parallel to 
the x-axis, hence 

x = x . 

i i 

y/ = y po + (yi _ y P o) cos ( z i _ z p 0 ) sin ^ 

z i' = z PO + (y i -y po ) sin<1> + ( z i - z po ) c °s<d 

where is the angle of rotation. (See Figure 3). 

Introducing the Notation: 

Ax. = x A -x po with corresponding definitions for Ay. , Az if equation 
(4) can be rewritten: 

*i' = x i 

Ay = Ay . cos 4> - Az . sin $ (5) 

A z / = Ay i s in 0 + Az . cos $ 

Equations (4) or (5) give the coordinates of any point attached to the 
paddle in its rotated position. 

The rotated position coordinates of any three points on the 
surface can be used to obtain the vectors describing the planes 
containing the paddle faces. It is, however, more convenient to 
carry out the transformation directly on the vectors , B., C i ,D.J 
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describing the plane, by equation (2). The 4-vector • C. , Dj ) 

describing the rotated plane is given in terms of Xj^, Yj^, Z{ mn by 
means of equation (3) where 

with corresponding definitions for Yj' mn and Z l ' mn . 

Further 

AX = X - x 0 

po 



AY = Y - y U 

po 

(6) 


AZ = Z - z 0 

po 


where the subscripts l. 

m,n have been dropped for convenience. 

In 

terms of these quantitie 

s: 


X 

= X 


Y' 

= y U + AY cos <p - AZ s in * 

J po 

(7) 

Z' 

= z U+AYsin<I> + AZcos<l> 


po 


Substitution of Equation (7) into (3) then results in the following equation 

A' = 

A 


B' = 

B cos <t> - C sin $ 


C' = 

B sin $ + C cos <5 

(8) 

D' = 

D - B [y po (cosd>- l) + z po sin*] 



+ C [ y p° sin$_z P o (cos* - 1)] 
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3. OPEP Points. 


The OPEP is assumed to be carrying only experiments, hence only 
the formulae concerning transformation of points are necessary. Further, 
the rotation takes place about an axis parallel to the z-axis. The formulae 
are consequently: 


X.' 

1 

" X To + 

( X i “■ X To) COS * " 

<< 

h*- 

1 

<< 

H 

0 

sin m 


yi' 

= y To + 

( X i " X To) sin,P + 

1 

<< 

0 

cos 

(9) 

z 

i 

- z . 

i 






where ¥ is the angle of rotation. (See Figure 4). 

4. Sunny or Shady Side Test. 

In terms of these vectors the sign of the product, 

• R^ (s. • R k j = s (See Appendix B) (10) 

determines the position of points r. and r k relative to surface S { . In 
particular if s< 0, r £ and r k lie on opposite sides of the surface, and if 
s > 0, r . and r, lie on the same side of the surface and for s = 0 at 
least one of the points lies on the surface. If now r . represents the 
Sun's coordinates and r k the coordinates of an experiment to be tested, 
s > 0 means then that the surface S i cannot cast a shadow on experi- 
ment r k . If, however, s < 0 further tests have to be performed to 
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determine whether the experiment lies inside or outside the shadow 
cylinder of the face. 

For the box surface S x is defined by equation (3). For a solar 
array surface Sj is defined by the 4-vector (A', B', C',D') which is 
given in equation (8). 

If an experiment is attached to the main box its augmented position 
vector 


R i = (*i • yj . *j • (ii) 

may be substituted directly into equation (10) for the test. If the ex- 
periment is attached to the solar array 

5 i ■ (“i' • yj’- */• ») U2) 

where 

x/.y/, z; 

are given by equation (4). If an experiment is attached to the OPEP 

gj ■ */• i) (13) 

where 

x/.y/, 2j ' 

are given by equation (9). 
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5. Shadow Cylinder Test 


If s < 0 (where s is defined in equation (10) above) this establishes 
that the Sun and the experiment lie on opposite sides of the surface. 
Further tests have to be performed to determine whether the experiment 
lies inside or outside the shadow cylinder of the face. 

This testing is accomplished by projecting both face and experi- 
ment onto a plane perpendicular to the Sun's radius vector. A polygonal 
face will always be projected into a polygon, unless the Sun's radius 
vector is parallel to the face, in which case no further testing need be 
performed for this combination of surface and experiment. Le t the first 
vertex r ^ (chosen arbitrarily from among the M vertices on the face) and 

the other vertices r p , r q etc. (taken in order around the face) 

be renamed r i (i = 0 . . . . M - l). 

If r o is now chosen as origin, the coordinates of the vertices 
relative to this origin are given by: 

*7 = 7 i - 7 (1 4 ) 


in particular 


r o = ( 0 , 0 , 0 ) 


Further 


r 


sv 



(15) 
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where r sv is the vector from satellite to Sun. (Note: because of the 
relative size of the components of ? sv and r o the vectors F s ' v and r gv 
will in general be numerically identical). 

In terms of these vectors, a coordinate system is now set up in 
the test plane (i.e. the plane perpendicular to 7 S ' V by means of the 
unit vectors i", j " . 


r sv x r i 
j r sv x r j 

/ • H 

r sv x 1 


(16) 


The coordinates of the projections of the vertices onto the test plane 
are now given by: 


p; - 


Yi J 


(17) 


where 


x. 

i 


■ i/ 

* l 



(i = 0 ••• M-l) 
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and the coordinates of the projection of the experiment (n) are given by 


q " = x " i " + y " i 

“ n n * n J 


// — . / • n 

“ r n • 1 

n n 


r n ' J 


In this test plane, the shadow cylinder is represented by a polygon with 
the M sides 


L jlj+ x (j = 0 •••M-l) 


where L j . +1 is the straight line containing the projection of the 
vertices r / and r j' +1 (Note: In this connection it is convenient to 

equate r ra ' to r o '). The equation of the side j+1 is given by (See 
Appendix C): 


or alternately: 


tj x* + b. y" + Cj = 0 
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where 


a i - yj 


j+i 


h. = - x " + x " 

Dj X j T X j+ , 


C = X " V " - X " V " 

c j x j yj+i x j+i yj 


( 21 ) 


The vector b^, c ^ is taken to represent the line j+I and the 

three vector Q n = (x n ", y n ", l) represents the point q n "; and let 


P 


j + 2 


(« 


II 

j+2. 



represent the j + 2 vertex of polygon projected onto the test plane. In 
terms of these vectors, the experiment represented by the vector Q n 
lies inside the polygon in the test plane (and consequently inside the 
shadow cylinder of the face) if: 


u = (Q. ' Lj.in) ( p i.2 L i.j.i) > 0 (22) 

for all j between 0 and M- 1. See Appendix D for the mathematical 
derivation of this expression. 

C. Ideal Control Laws and Definition of Axes Systems 

The geometry of the OGO satellite is described in a coordinate 
system attached to the main structure, (box coordinates), whereas the 
orbit and consequently the vectors from the Sun and Earth to the sat- 
ellite as well as the velocity of the satellite are given in a different 
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coordinate system. In order to carry out the test described in section 
IVB, the components of all vectors have to be known in a single coordinate 
system. It is convenient to use the box coordinate system to carry out 
this testing, since fewer transformations have to be made. 

In order to relate these coordinate systems and to establish the 
transformation matrices, the ideal control laws for the OGO satellite 
were assumed. Further, these relations are written entirely in terms 
of the vectors r VE , 7 yE and r gy . Consequently, the transformations will 
be valid, no matter what coordinate system is used in describing the 
orbit. 

In order to describe the transformation from the orbital to the 
box coordinate system, a matrix is established. In addition, an angle 
of rotation must be determined for the solar array and another for the 
OPEP. 

In the following, four coordinate systems are considered. Each 
coordinate system is defined by a triple (i, j , k) of mutually orthogonal 
unit vectors. The systems attached to box, solar array and OPEP are 
identified by subscripts B, p and T respectively. Quantities without 
subscripts refer to the orbital coordinate system. 

The three vectors r yE , r yE , r sy (the position and velocity vector of 
the satellite with respect to the Earth and the position of the Sun with 
respect to the satellite) are given in the orbital coordinate system. 


27 



1. The Orbital To Box Transformation Matrix. 


The ideal control laws are specified in Reference 5. 

a) . The positive z B face is aligned with the radius vector from the 

center of the Earth to the satellite. 

b) . The solar array is normal to incident sunlight. 

c) . The Sun does not shine on the positive y B face. (Reference 5). 

d) . The positive x face of the OPEP faces forward in the plane of 

the orbit. (Reference 6). These laws lead to the following 
equations : 

Due to a). 


k 


B 



Due to b). 

r sv 
r sv 



(23) 


(24) 


Since the paddle rotates about an axis parallel to i B 

j p = j B cos$» + k B sin$ (25) 

From (equation 25) 


jp 



0 
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and by definition 


k 


B 



= 0 


Consequently 



The vector triple is completed by defining 

1b ~ x *B 


According to c). 


1b 


sv 


< 0 


Equation (28) resolves the ambiguity in equation (26) as follows: 
stituting from (26) into (27): 



From (24) 


r 


sv 


ip 



(26) 


(27) 


(28) 

Sub- 


(29) 


(30) 


29 



From equations (28), (29) and (30) 


0 > i B ' r sv 1 | r sv|[(^B‘ jp) 2 k B 2 jp 2 ] (31) 

From (25) k B • j p = sin$ 

0 > T b * r sv = ± |7 sy |[sin 2 $-l] (32) 


In Equation (32) the term |7 SV | [sin 2 $>- l] is less than or equal to zero, 
it follows that the plus sign must be chosen in equations (26) and (29). 

If now r is an arbitrary vector with components x, y, z in the orbit 
system and coordinates x B , y B , z B in the box system 


x 


B 



y B 




(33) 


or in matrix notation 



(34) 


The vectors i B , j B , k B are defined by equations (23), (26) and (27). 
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2. The Solar Array Angle $. 


The solar array angle $ is obtained from equations (24) and (25). 

r sv — 

• j B = cos$ 

(35) 

• k B = sin<I> 

For the transformations, only cos 4> and sin$ are needed, the angle 
need not be computed. 

3. The OPEP Angle <F. 

The vector i T is in the plane of the orbit and 



k 


T 



(36) 


The orbit plane is defined by the angular momentum' vector H 


H r VE x r yE 


(37) 


H is normal to the orbit plane hence 

_ k T x H 

i_ = ± Tn rrr 

| k T xH | 

Further, i T points in the direction of motion, i.e; 


(38) 


( 39 ) 
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Substituting (23) and (37) into (38) yields: 


“ ± 


= ± 


( r VE) x ( r VE x r VE ) 


r yg k T x H 


7 k x h| FyE * ryE ) FyE ^ ryE " ryE ^ ryE l 

VE T 


( 40 ) 


Substituting (40) into (39): 


0 < 


^ 1 

r,- 1 2 I'i: i 

JL * r r> ± - — 

T VE r VE k X xH 

y r VE 1 I r VE 1 


by Schwartz's inequality the expression in braces is positive, hence 
the plus sign has to be selected in equation (38). 

The vector triple is now completed by: 


j f x i x 

The OPEP angle ¥ is now obtained from the equations: 


(41) 


i T • i B = cos? 

i T ' T b = si"'* 1 


(42) 


Again, for the transformations, only cos'? and sin Ware needed. 

D. Detailed Shadow Testing Procedure. 

The equations of the planes containing the potential shadowing 
faces are established initially and the corresponding vectors (Aj, 
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B it Cj, Dj,) . As described in equations (1) and (2) are computed and 
stored, together with the list of its vertices. This is done initially for 
both box and solar array faces. The solar array faces, however, must 
be rotated at every test time after the angle 4> has been established 
according to equation (8). The testing at every point is now carried out 
in the following sequence: 

1. Satellite Earth Shadow Testing. 

It is determined whether the satellite is in the Earth's shadow 
cone. If it is, the possible shadowing effect of satellite components 
need, of course, not be determined. If this test, however, shows the 
satellite to be in the sunlight, testing has to be carried out for each 
individual experiment as follows: 

2. Individual Experiment Testing. 

a) . The matrix described in equation (34) is computed. 

b) . The solar array angle 0 is computed by equations (35). 

c) . The OPEP angle 'P is computed by equations (42). 

d) . The coordinates of Sun and Earth with respect to the 

box are transformed from the orbital to the box 
coordinate system by equation (34). 

e) . The coordinates of the experiment are transformed, 

if necessary, (i.e. if the experiment moves with either 
the solar array or the OPEP) by equations (4) or (9) 


respectively. 
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Items a) through d) of this list to be performed only once for each set 
of tests (i.e. for each time point). Item e) has to be performed for each 
experiment. 

After all this setup work has been completed, the shadow tests are 
carried out for this particular experiment designated by (n) as follows: 

All the faces are tested in turn until: 

a. Either all surfaces have been tested and it has been established 
that none casts a shadow on experiment (n) . In this case, experiment (n) 
has been established to receive incident sunlight and the next experi- 
ment (n+ 1) is tested. 

b. Or a face has been found which casts a shadow on experiment 
(n) . In this case, experiment (n) has been established not to receive 
incident sunlight, the effect of the remaining faces need not be examined, 
and the next experiment (n + 1) is then tested. 

The hierarchy of tests, for a particular experiment (n) and face (i) 
is described in the following: 

1. Test whether face (i) is connected to solar array, and carry 
out, if necessary, the transformation described in equations (8). 

2. Test whether the experiment lies on the sunny side or the 
shady side of the plane containing face (i) by using the augmented 
box coordinates of the Sun as IT and the augmented box coordinates 
of the experiment as R k and substituting in equation (10). 
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If s > 0 the experiment is on the sunny side, this face cannot cast 
a shadow on experiment (n) and the next face is then tested. 

If s = 0 either the Sun, or the experiment, or both lie on the sur- 
face and the surface is assumed not to cast a shadow on experiment (n) . 
In this case, the next face is tested. 

If s < 0 the experiment (n) lies on the shady side of face S. . See 
Appendix B for more details. If it also lies within the shadow cylinder 
of the face it will be in shadow. 

3. Test whether experiment (n) is within the shadow cylinder of the 
face (i). In order to carry out this test, the following steps have to be 
taken: 

a) . Compute the vectors r.', r ' and r s ' v for the experiment 

(n) , all the vertices of faces (i) and the Sun by equa- 
tions (14) and (15). 

b) . Set up the coordinate system ( i", j") in the test plane 

by equation (16). 

c) . Compute the coordinates (x", y" ) of the vertices and 

the experiment in the test plane (i.e. the plane normal 
to the Sun's radius vector), by equations (17) and (18). 

d) . Compute u by equation (22) for edge (j). 
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If u < 0, the experiment lies outside the shadow cylinder of this 
face, this face does not shadow experiment (n) and the next face i + 1) 
is then tested. 

If u = 0 the experiment is also taken to be outside the shadow 
cylinder and the next face (i + 1) is then tested. 

If u > 0 for the j th edge the next edge (j + 1) of the face is tested 
until one edge has been found for which u<0. Then this face (i) 
does not cast a shadow and the next face (i + l)is tested. If all 
edges of a face have been tested and u> 0 for all edges then this 
face (i) casts a shadow on experiment (n) . The program then 
proceeds to the next experiment (n + 1). This is explained in more 
detail in Appendix D. 
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V. MATHEMATICAL MODEL FOR HEAT INPUTS TO AN 


ARBITRARILY ORIENTED SURFACE OF OGO. 

A. Direct Solar Radiant Heat. 

This is the heat flux q ns (t) incident upon a surface of the nth ex- 
periment at time t, which radiates directly from the Sun. It is given by: 


q„ s (t) = I(t)af(i n (t))[|r sv (t) 


(43) 


here: I(t) = 0 or 1 depending on whether the experiment is in shadow 

or Sun, cris the incident radiation on a unit area at a distance of 1 AU 
from the Sun, and i n the angle of incidence given by 


cos i 

n 



(44) 


where n n is the normal to the experiment surface in box coordinates. 
If cos i n < 0 , the experiment surface faces away from the Sun and 
receives no radiation. This is expressed formally by 

*(*„) = l( COsi n + | COsi n|) (45) 

In the case of an earth satellite the variation of j r sv | is too small to 
matter and a constant value may be used. The accumulated radiation 
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between times t D and (t Q + N&t) is given by: 


% s (to- t o +NAt) - £ \ ["Ins (to +kAt ) (t Q + (k+ l)At)] A n At k (46) 

k= 0 

where 

A n = the area of experiment (n) . 

B. Reflected Solar Radiant Heat. 

This is the heat flux q na (t), incident upon a surface of the nth ex- 
periment at time t, which originates at the Sun and reflects from the 
Earth. In integral form (See Figure 5). 


%a(t> 



(47) 


where, 

a = the albedo = the fraction of the solar constant a which is 

reflected away from the Earth. The value for a is about 0.34 
(Reference 7). 

a = 1.97 ± 0.01 cal cm -2 minutes -1 = 0.1374 watts cm -2 = flux 
of total radiation received outside the Earth's atmosphere 
per unit area at mean Sun-Earth distance (Reference 7). 
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f a (^l > > ^3 » F ) 


cos y x cos y 2 cos y 3 


77r" 


if cos y l > 0 and 
cos y 2 > 0 and cosy 3 >0 


f a = 0 if cosy x <0 or cosy 2 <0 or cosy 3 <0 
dA = an elemental area of the Earth's surface (Figure 5). 
y 1 = the angle of incidence of the Sun's rays on the Earth's sur- 
face element (Figure 5). 

y 2 = the angle between the radius vector r which points from dA to 
the vehicle and the normal to the Earth's surface element dA 
(Figure 5). 

y 3 = the angle between the vector r and the unit outward normal 
to the experiment surface (Figure 5). 
r = the vector which points from the vehicle to dA (Figure 5). 

The method of numerically evaluating Equation (47) is given in the 
following paragraphs. 

Let 6 be the angle between r yE and the radius vector from the 
center of the Earth to the surface element in question, and a be the 
angle corresponding to the longitude in the plane normal to r yE the 
zero line being the projection of i B on this plane. (See Figure 5). 

The part of the sphere which the satellite can see is then 
characterized by 

0 ! a < 2rr 

0 <6 < jB 
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where 


a E 

C ° S /3 = T~ j (48) 

I 1 VE I 

here a £ is the radius of the Earth. The range of 8 and a is now sub- 
divided such that 

A0 = yS/L 
Aa - 2 tt/M 


The normal to the Earth 1 s surface element AA, is approximated by 


N U “ i B sin^j cosa m - jg sin(9j Sina m - k B COS0J (49) 


where 


(i + JV 


l = 0, 1, 2, 3 ’■* , L-l 


2n 


(■*i) (50) 

’ m - 0, 1, 2, 3 , M- 1 


M 


The position vector of this surface element from the satellite is given 
by: 


r 


VE 


a E N in 


(51) 
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Note that the r VE vector in this equation is given in the box coordinate 
system. 

It may now be determined whether radiation from this segment of 
the Earth's surface reaches experiment (n) by using vector r lm in place 
of the Sun vector r gv in the tests described in section IV. Note in 
particular that the only additional set up work which has to be done is 
the transformation of the experiment into box coordinates (item e when 
necessary) and only the tests have to be performed. 

The cosines of three angles determine the amount of incident 
radiation received from the lit surface of the earth. They are: 

a). The angle of incidence y 1 of the Sun's rays on the surface 
element. This cosine determines the intensity of the 
illumination of the surface element and is given by the 
equation. 


cos ?lU 




( 52 ) 


where 


S 1 01 


r 


SE 


a E ^ lm 
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The degree of accuracy required for this computation allows 
7 , = 7 CS . and equation (52) then becomes: 

— r SE 

co s y i j m - N lm • (53) 

| 1 SE | 

b). The angle y 2l between the radius vector from the sur- 
face element to the vehicle and the normal to the sur- 
face element. This angle controls the size of the surface 
element as seen from the satellite. The cosine of this 
angle is given by equation (54). 


cos y 


2 lm 



N 


lm 


(54) 


r lm in this equation is defined by equation (51); 
c). Finally the angle y 3lm between the incident radiation and 
the normal to the experiment surface determines the solid 
angle subtended by the experiment at the origin of radia- 
tion. This angle y 3 j m is given by equation (55) 


cos y 


3 lm 



(55) 


n n here is the normal to the experiment surface in box 
coordinates. 
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1. Contribution due to a lit surface element: 


The contribution of the particular element AA^ m to the reflected 
solar heat flux at time t is given by equation (56). 


^„a(0 


aaI tn, f (?2 J ( f 3 Jm) 


77 r 


T I 2 


Im 


( 56 ) 


T 0 if the experiment is shaded by the satellite 
I 1 if the experiment is not shaded 


Here 


f(y) 


= \ (• 


cos y + I cos y 


l) 


and 


AA im - a £ 2 A0 Aa sin 0 ^ (57) 

In equation (56)7 Un)t 7 2Jm , 7 3lm and r are functions of time. The 
total reflected solar heat flux received from all surface elements is 


^na^) = ^"a (t) 

l , m 


(58) 
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This may be accumulated over a time interval from t Q to (t Q + NAt) 


^ | 

= ( t o + < k+ 1 > At )] A n At k (59) 

k = 0 

C. Earth Emitted Radiant Heat 

This is the heat flux q nE (t) incident upon a surface of the nth ex- 
periment at time t, which is emitted by the Earth. In integral form 
(See Figure 5). 


*nE 


(t) 


A 



y 3 , r) dA 


(60) 


where 


A = 0.02215 watts cm” i 2 = the emittence per unit 

area of the surface of the 
Earth. 


cos y 2 cos y 3 


77T 


i f cos y 2 > 0 and cos y 3 > 0 


= 0 


i f cos y 2 = 0 or cos y 3 = 0 


This is evaluated numerically and is analogous to (56), (58), and 
(59), except that the angle of incidence of solar energy is immaterial. 
The equations are as follows: 
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The contribution from the element AA, to the Earth emitted heat 

i m 

flux at time t is given by equation ( 61 ). 


^nE (O 




The total Earth emitted heat flux received from all surface elements is 


(t) 


2 >„, 


The accumulation over a time interval from t o to (t Q + NAt) is 


n- i 

(«.■ *. tNit ) = 2l 1 [ 4 " e (*• +kAt ) (*• * (k + 1)At )] K At k (63) 
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VI. RESULTS AND DISCUSSION 


A. Shadow Study 

Figures 6 through 8 present the amount of time per orbit that the 
EGO S-49 experiments spend in the shadow of the earth, the main 
satellite structure (the box) or of the solar array. These data are for 
three faces of the main box +z B , -z B and -y B for the boom or appendage 
mounted experiments (E. P. 1 through E. P. 6), and for the orbital 
plane experiment package (OPEP 1 and 2). The coordinates for the 
sides of the main box and the experiments are shown in Section III of 
this document. In the case of the shadowing of the faces of main box 
test points were assumed to be along the axes of the box coordinate 
system and slightly above the surface (0. 1 inch). The shadow history 
of the test point of each face is assumed to closely approximate the 
shadow history of the face itself. The test point device was incorporated 
to facilitate the mathematical description and programming of the 
system. The coordinates of the main box test points and normals are 
as follows: 
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Coordinates 



X 

y 

Z 

Box Experiment 

(inches) 

(inches) 

(inches ) 

+x 

14. 707 

0 

0 

-x 

-14. 707 

0 

0 

+y 

0 

23. 6 

0 

-y 

0 

-43. 6 

0 

+z 

0 

0 

15. 083 

-z 

0 

0 

-15. 083 

Normals 




+x 

1. 0 

0 

0 

-x 

-1. 0 

0 

0 

+y 

0 

1. 0 

0 

-y 

0 

-1.0 

0 

+z 

0 

0 

1. 0 

“Z 

0 

0 

-1. 0 


The study for shadowing is based on a nominal orbit (See Section III for 
injection conditions) with a launch date of May 15, 1964 at 13. 0 hours 

U. T. 
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1. The Main Box Experiments 


Figure 6 shows the history of shadow time per orbit for a year of 
flight time for the +z g , -z B and -y fi faces of the main box. The figure 
indicates that for the +z g and -z fi faces the shadow time per orbit 
varies between approximately 1 hour and 41. 0 hours. The -y fi face 
experiment, however, spends practically all of its time in sunlight. 
During the first 34 days in orbit the -y B face experiences shadows of 
about 0. 5 hours after which time it senses full sunlight except between 
the 144th and 152nd days after launch when the shadow time gets only 
as high as 2.5 hours. At about the 265th day after launch, the -y B face 
again goes into a period during which it senses about 0.5 hours shadow 
per orbit. The shadow history of this face (in fact of the whole satellite 
and experiments) is a function of the time of launch. Consequently, for 
other launch dates the shadow history will be different. 

2. Boom Mounted Experiments 
The time -in -shadow history of the boom mounted experiments 
(E. P. 1 through E. P. 6) is shown in Figure 7. This time history curve 
is typical of any launch date, however the curve set would be different 
for different launch dates. 

The figure clearly shows that the experiments in the order E. P. 1, 
E. P. 2, E. P. 3, E.P. 5, E. P. 6 (E. P. 4 = E. P. 6) are exposed to 
increasingly more sunlight. The maximum shadow time per orbit 
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varies from almost 40 hours for E. P. 1 to about 2. 5 hours for E. P. 4 


and 6. 

3. OPEP 

Figure 8 presents the time -in-shadow histories of the orbit plane 
experiment packages (OPEP). The experiment packages, OPEP 1 
and OPEP 2 experience some shadow time during every orbit. These 
shadow times, however, vary from 41. 0 hours per orbit to less than 
1. 0 hours per orbit. Again, as in the other cases, this curve set will 
be different for different launch dates. 

B. Heat Inputs 

Figures 9 through 20 present the heat inputs incident upon the 
different faces of the satellite as a function of time from perigee for 
the first orbital period (42. 75 hours). These heat inputs include the 
Earth emitted heat flux q £ , the direct solar heat flux q s and the solar 
reflected heat flux q a . The solar heat input to an experiment is set 
to zero if that particular experiment is shielded from the Sun by the 
Earth, the main box or the solar array. If an experiment is in the 
eclipse of the Earth the reflected solar (as well as the direct solar) 
heat input is set equal to zero. However, the effects of shadowing of 
an experiment from the Earth by the main box or the solar array is not 
considered in the computation of the solar reflected or Earth emitted 
heat inputs to that particular experiment. The reflected solar heat 
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flux is calculated by two different methods. The exact value is com- 
puted by numerically evaluating the integral given by equation (47). 

The method of calculating the approximate values is given in Reference 
8. In this method an average value is assumed for y l in equation (47). 
It is assumed that y 1 (average) is the angle between the vector from 
the Earth to the Sun and the vector from the Earth to the satellite. In 
other words for a given location of the Earth, Sun and vehicle, the 
approximate reflected solar heat input depends only upon the angle 
between the experiment's normal and the Earth vehicle vector. The 
Earth emitted and solar reflected heat inputs are negligibly small 
when the satellite is beyond an altitude of 20, 000 nautical miles or 
correspondingly 2. 63 hours of flight time from perigee. The data in 
these figures are based on a time of injection into orbit of May 15, 

1964 at 8. 0 hours U. T. These curves are intended to demonstrate the 
capability to calculate heat inputs rather than to provide design data. 
The values of the Earth emitted heat inputs depend only upon the shape 
of the orbit and are independent of launch time. The direct solar and 
reflected solar heat inputs depend upon the relative locations of the 
Earth, the Sun and the vehicle; consequently they are dependent upon 
the launch time. The unit normals to the various S-49 experiments 
are given in section III. 
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1 . The Main Box 


a. Earth Emitted Heat Inputs 

Figure 9 presents the Earth emitted heat inputs q E to the main 
box as a function of time from perigee. The + z b face has the largest 
Earth emitted heat input since it always faces the Earth. Conversely, 
the — z b face has a zero Earth emitted heat input. The +x B , -x Q , +y B , 
and the -y b faces all have equal Earth emitted heat inputs due to geo- 
metrical symmetry with respect to the Earth. 

b. Direct Solar Heat Inputs. 

Figure 10 presents the direct solar heat inputs (q s ) to the main 
box as a function of time from perigee. This is a typical curve. The 
solar heat inputs to the +x B , -x B and +y B face are always zero. This 
is true since the OGO satellite is controlled such that these faces 
receive no direct solar heat. 

c. Reflected Solar Heat Inputs. 

Figure 11 presents a comparison between the exact and approxi- 
mate values of the reflected solar heat inputs. The +z B face has the 
largest reflected solar heat input since it faces the Earth. The -z B 
face has zero reflected solar heat inputs since it faces away from the 
Earth. The exact value of the reflected solar heat inputs to the +x B 
and -x B faces are equal since there is symmetry with respect to these 
faces and the sun-lit portion of the Earth. The approximate values of 
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the reflected solar heat inputs to the +x B> -x B , +y B and ~y B faces are 
all equal. This is true since an average value is assumed for y x in 
equation 47. It is assumed that y (average) is the angle between the 
vector from the Earth to the Sun and the vector from the Earth to the 
satellite. This is tantamount to assuming that for a given Earth, Sun 
and satellite configuration the amount of reflected solar heat inputs 
depends only upon the angle between the normal to the surface and the 
Earth -satellite line. Figure 11 shows that the largest difference 
between the exact and approximate value is 1. 16 milliwatts /cm 2 . 

2. Boom Mounted Experiments 

Figures 12 through 14 present the heat inputs to the boom mounted 
experiment packages (E. P. 1, E. P. 2, E. P. 3, E. P. 4, E. P. 5, and 
E. P. 6) as a function of time from perigee. The normals which define 
the orientation of the faces which have the various heats incident upon 
them are given in Section III of this report. These particular faces 
were chosen for investigation since they have the largest openings through 
which radiant heat may enter to the experiment packages, 
a. Earth Emitted Heat. 

Figure 12 shows the Earth emitted heat inputs to the experiment 
packages. All six experiment packages have equal Earth emitted heat 
input time histories. This is true since the normals to these experi- 
ments are parallel to the surface of the Earth. 
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b. Direct Solar Heat Inputs. 


Figure 13 shows the direct solar heat inputs to the experiment 
packages. E. P. 4 and E. P. 6 have the same value of solar heat inputs 
since their normals point in the same direction. The direct solar heat 
inputs to E. P. 1, E. P. 2 , E. P. 3 and E. P. 5 are identically zero. 

This is true since the normals to these experiments are oriented such 
that they never receive sunlight. 

c. Reflected Solar Heat Inputs. 

Figure 14 presents a comparison between the exact and approxi- 
mate values of the reflected solar heat inputs to the boom mounted 
experiment packages. E. P. 4 and E. P. 6 have equal exact reflected 
solar heat flux since their normals are equal. Likewise E. P. 2 , E. P. 3, 
and E. P. 5 have equal exact reflected solar heat flux because their 
normals are equal. E. P. 1 through E. P. 6 have equal values of 
approximate reflected solar heat flux since the angles between their 
normals and the Earth vehicle vector are all equal, i. e. , 



Figure 14 shows that the largest difference between the exact and 
approximate values is 1. 36 millwatts / cm 2 . 
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2. The Solar Oriented Experiment Packages 


The solar oriented experiment packages {SOEP 1 and SOEP 2) 
both have their normals pointing toward the Sun and therefore have 
equal heat inputs. 

a. Earth Emitted Heat Inputs. 

Figure 15 shows the Earth emitted heat flux to SOEP 1 and 2. 

b. Direct Solar Heat Inputs. 

The direct solar input to SOEP 1 and 2 is constant (137. 4 milliwatts/ 
cm 2 ) since it is controlled to point toward the Sun. Figure 16 shows 
that the satellite did not go into an eclipse during this particular time 
period, since the curve is continuous. 

c. Reflected Solar Heat Inputs. 

Figure 17 presents a comparison between the exact and approxi- 
mate values of the reflected solar heat inputs to SOEP 1 and 2. Figure 
17 shows that the largest difference between the exact and approximate 
values is 1. 22 milliwatts / cm 2 . 

3. The Orbit Plane Experiment Packages . 

The orbit plane experiment packages (OPEP 1 and OPEP 2) are 
controlled so that their normals lie in the orbit plane and such that 
the angle between their normals and the vehicle velocity vector is less 
than 90° (or n • r y£ > 0 ). 


54 



a. Earth Emitted Heat Inputs 


Figure 18 presents the Earth emitted heat inputs to OPEP 1 and 2. 
The values are equal since both OPEP normals point in the same direc- 
tion. 

b. Direct Solar Heat Inputs 

Figure 19 presents the direct solar heat inputs to OPEP 1 and 
OPEP 2. The two experiments have the same inputs except for shadow- 
ing of the OPEPs by the main box. This figure shows that OPEP 2 has 
no solar heat input from 0. to 32. 5 hours from perigee. It also indi- 
cates that OPEP 2 is shadowed from the Sun by the main box. OPEP 1 
becomes shadowed by the main box at 42. 5 hours. 

c. Reflected Solar Heat Inputs 

Figure 20 presents a comparison between the exact and approxi- 
mate values of the reflected solar heat inputs to OPEP 1 and 2. OPEP 
1 and 2 have equal inputs since they have equal normals. The largest 
difference between the exact and approximate values is 0. 7 3 milliwatts/ 
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APPENDIX A 


The Equation of a Plane 

Equations (1) and (2) in Section IV -B show a fourth order determinant 
and the analytical expression of a plane surface. The determinant and its 


relation to the analytical expression are shown here. 


Letrj , r B , r n be vectors in an orthogonal coordinate system which 
are directed to any three vertices of a plane surface. 

If i, j, k are unit vectors along the x, y, z axes respectively r ; , 
t m , r n can be expressed as: 


i-j - ixj + jy j + kz l 

r = ix + jy + Icz 

m m m m 

7 = lx + jy + Icz 

n n ■'•'n n 

r = ix + jy + kz 


(A-l) 


Where r is any vector from the origin to the plane surface. Vector 
subtraction forms the vectors u, v , w in the plane: 


u = 7-7j = (x-x l )i + (y-y I )j+(z-z l )k 

V = 7 - 7 m a (x -X m ) i + (y- y m ) J + (z - zj k (A-2) 

w = 7 - 7 n = (x-x n ) i + (y-y n ) T + (z-z n )k 
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Since u, v, w are coplanar their triple scalar product is equal to 


(uxv) • w 


u 

u 
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X 
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V 

V 

V 

X 
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w 

w„ 

w 

X 

y 



z 


z 


z 


0 


Where the elements of the determinant are the coefficients of i , 
in equation (A- 2). Hence 


( x-x z) 

(y 

-».) 

( z ' z l) 

( X - X m) 

(y 

-y.) 

( Z - Z m) 

( X " X n) 

(y 

-o 

( Z " Z n) 


The above determinant may be expanded to: 
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zero. 


(A-3) 


j , k 


(A-4) 


(A-5) 
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Factoring out -1 from the l, m, n rows of each of the above determinants 


and interchanging the rows for symmetry yields: 


x y z 


x y z 

x y z 

m J m m 

- 

x i y i 

x y z 


X y z 
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(A-6) 


The determinants of equation (A-6) can be considered the minors of an 
expanded fourth order determinant which has been expanded for a 
column whose elements are unity. 



(A-7) 


Hence, the expression of a plane is described as a fourth order 
determinant. 
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Expanding equation (A -7) for the elements of the first row produces: 
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Letting the minors of x, y, z, 1 equal A, B, C, D respectively the 
analytical expression of a plane surface can be written. 


(A-8) 


Ax + By + Cz + D = 0 


Accordingly, any plane surface can be described by the proper choice 
of A, B, C, D. For instance the i th surface is: 


A. x + B. y + C. z + D. - 0 (A-9) 

or in other words the four vector, (A. , B , C , D.) describes the 1 th 
surface. Since a triple scalar product can be written as a third order 
determinant, the coefficients of the analytical expression for a plane 
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surface can be written as: 


A. 

1 


B. 


C 


where 


y i 

y m 

y„ 


l 

l 

l 

l 

l 

l 






= U-(X u „xY,.J 


y„ 

y„ 


^Imn (^Jmn X ^Jmn) 


X Un = ( X i- X m- X n) 

Y i mn s (y i • y m > y n ) 

Z lmn S ( Z l> Z m’ Z n) 


u - (1, 1, 1) 


(A-10) 


(A-ll) 


The right hand sides of these definitions are three-vectors. 
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APPENDIX B 


The sign of a mathematical expression may be used to determine 
whether any two points are on the same or opposite sides of a plane 
surface. 


where 


if 


(Si • R^ (Sj • R k ) = s Equation (10), Section IV 


S. 


R. 


is a four -vector denoting the 
analytic representation of a 
plane containing the i th face. 
(See Appendix A). 

(xj , yj , z. , 1) is a four -vector denoting the 

position of the Sun. 

( X k i y k . z k , l) is a four -vector denoting the 

position of an experiment 
package. 


s >0then R. and R. are on the same side of S.. 

J K 1 

s <0then R. and R k are on opposite sides of S. . 
s = Othen either R. or R k , or both are in the plane of S. . 


In general, to determine the relative positions of two points in 
space with respect to a surface it is sufficient to establish the angles 


B-l 


between a normal vector to the surface and the position vectors 
measured from the surface to the points. 

Graphically this may be represented (See Figure B-l), where 

n = normal vector 
q j = vector from the surface to r . 
q k = vector from the surface to r k 
r = position vector of the Sun 
= Xj I + y i j + Z. It 

r k = position vector of an experiment package 

= x k 1 + Vk T + z k k 
r = vector to origin of n 

= xi + yj + zk 

<& = angle between n and q. 

<p k = angle between n and q k 

Note that n, q. , q k , r . , r k and r are all three -vectors (i. e. defined in 
3-dimensional orthogonal space), and that 


r i t R , 
r k / R k 
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Based on these definitions the following relationships are made: 


CASE 


A 

SIGN_OF 
n • q. 

SIGN OF 
n *q k 

SIGN OF 

(n-qjx 

( n ‘ %) 

RELATION OF 
q. andq k 

1 

<t> } <rr/2 

4> k < v/2 

+ 

+ 

+ 

same side 

2 

4>. >tt/2 

<£ k <rr/2 

- 

+ 

- 

opposite sides 

3 

4>. <n/2 

<t> k > tt/2 

+ 

- 


opposite sides 

4 

i 

4>. >tt/2 

CM 

> 

A 

- 


+ 

same side 


This table may be summarized in one mathematical expression, 

(n-i.)(n-q k ) = s (B-l) 

When s>0, q^ and q k are on the same side of the surface. 

When s < 0, q. and q k are on opposite sides of the surface. 

In order to prove that the sign of • R • j ( l S i • R k j yields the 
same information as • q^ • q k ^ it is necessary to show the 
relationship between the two expressions. 

A. x + B. y + C. z + D. = 0 (equation of plane surface) 

V(A. x + B. y + C. z+D.) = normal to plane surface (B-2) 

V (A a x + B a y + C z + D £ ) = A i i + B i j + C A k 
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Hence, 


n = Aj i + B. j + C k 


From the figure above it can be seen that 


q i = r i " r 


% = r k " r 


These substituted into equation (B-l) yield: 

[n • (7j -7)][n • (7 k -7)] = s 

[(n • - (n • 7)] [(7 • r k ) - (n • 7)] = s 

[( A i x j +B i Vj +c i z j) - ( A i X + B i y + c i z )] x 

[( A i x k +B i y k +c i z k) - ( A i x +B i y + c i z )] 

Since 


- s 


-D i - Aj x + B i y + C A z 


Hence 


( A i*i +B iy, *C l!j tD 1 )(»,« l tB l , 1 .C^tD,) = S (B-3) 
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Defining 


R i “ (wv 1 ) 

\ s (* k .y k . * k . *) 

and noting that 

s i * <*, . b , . c , . »i) 

Equation (B-3) can be written as 

(s, • *,) (s, • S k ) = . 

Hence, it is shown that the sign of this expression (which is identical 
to equation (10), Section IV) is a valid means of determining whether 
any two points are on the same or opposite sides of a plane surface. 
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APPENDIX C 


The Equation of a Straight Line 

The projections of the vertices of a polygonal face onto the test 
plane (See Figure C-l) are given by: 


ii 


= x . 


y ; j 


(C-l) 


The line connecting two successive vertices, p/' and p.'^ can be 
represented by the vector difference of these two vertices. 


‘i.i*! = Pi" - P/"*i (C-2) 

T /.i*. = h"-*A.) T " * (Vj-- vAi) T” (C-3) 

Since the above expression is the equation of a straight line of the form 

ax" + by" + c = 0 


where -a/b is the slope of the line. Hence, the slope of i. . + J is: 


Slope 


(A-yA.) 

( x j"- x Ai) 


(C-4) 


The slope of any straight line is the same anywhere along the line, 
consequently: 



j 


I*ilL 


- X 


j + 1 



(C-5) 
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VECTOR 
PERPENDICULAR 
TO POLYGONAL PLANE 


Rsv 




Expanding yields: 

-(*/ - v. .)(*'-»/) + 

-x."y" + x j " y + x." + 1 y" - x." + 1 y." + x" y." - x" y ." +1 


= 0 


c." y." + x." y" , 

i i j 3 j + i 


The above equation can be written in terms of determinants: 
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The determinants of equation C-6 can be considered the minors of an 
expanded third order determinant which has been expanded for a column 
whose elements are unity. 


x* 
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x j" 

yj" 

, it 

it 

i + 1 

y i+l 


= 0 
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Expanding equation C-7 for the elements of the first row produces: 

*•(*/- y.*.) - y(v-*A.) ♦ =° <c-«> 


Letting 


a . 
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b 
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“X 


V - y A. 

' + X * 


j + 1 

x." y " , - x." y." 

i j * l j + l J j 
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Hence, the familiar straight line formula: 


a j x " + b j y" + c j = 0 ' (C-9) 

This equation can be considered the result of a dot product of two three - 
vectors: 


L i , i + 1 “ ( a i ’ b j ’ c i) 

(C-10) 

II 

X 

H-* 

(C-ll) 

L j , i + i ‘ P “ aj x" + bj y" + Cj - 0 

(C-12) 
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APPENDIX D 


The Shadow Test 

In order to test whether or not an experiment is in the shadow of 
the satellite the experiment is investigated with respect to each surface 
of the satellite and to the cylindrical projected shadows from each of 
these surfaces. The polygonal surface and the experiments are pro- 
jected onto a test plane which is perpendicular to the Sun's light rays. 
(See Figure C-l, p. C-2. ) If a given point on the test plane (i. e. the 
projection of an experiment onto the test plane lies within the bounda- 
ries of the polygonal shadow on that plane the experiment will neces- 
sarily have to be in the cylindrical shadow since the Sun and experi- 
ment are on opposite sides of the satellite surface. If the boundary of 
the polygonal shadow is described by a series of vectors i. + J where 
j = 0, 1, 2, • • • M - 1 for an M sided polygon connected tip to tail in 
clockwise fashion and if the given point lies on the right hand side of 
each of these boundary —describing vectors the given point must lie 
within the boundary. If the point lies to the left of any side of the poly- 
gon the point cannot be within the boundary and hence is not in the 
shadow of that surface of the satellite. 
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Graphically this may be represented as: 



In this case the given point lies 
to the right of each vector and 
the experiment is in the shadow. 



Given Point 

In this case the given point lies 
to the right of some but not all 
of the boundary describing vectors 
and hence the experiment is not in 
the shadow of the surface causing 
this polygonal shadow. It may be 
in the shadow of another surface. 


The sides of the polygon are represented by: 


Z j,i+i ~ Pj"+ 1 “ P j" (See Appendix C) 

where 

j = 0, 1, 2, ••• M- 1 
M = number of sides on polygon 


and 


( 0 , 0 , 0 ,) 
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The vector q n " is the projection of an experiment package on the test 

plane in the p." system. 

Let 


E n - q n " ~ Pj" (See Figure C-l page C-2) 


and also let 

? ' Pj") xT i.i*l 

B = T*.,,,*? 

M lies in the plane of the polygon and is perpendicular to l. , + l . 

If M • E n >0 for all j from 0 to M - 1 the experiment projection will lie to 
the right hand side of each l . . + j and consequently in the polygonal 
shadow. If for any j , M • E <0 the experiment projection cannot lie 
within the polygonal shadow and consequently the experiment package 
cannot be in the cylindrical shadow projected from the polygon in 
space. It is shown here that 

M • E = (q • L , + 1 ) (p. . • L. . . 

n \ n j » j + 1/ \ j + 2 3 . j + 1 J 

Where 


’ kj j 


and P j+2 
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are three-vectors defined on page D-6 


ViM = = (V*.-*/) T ' + (y A.-yj") T” 

(vr ; ,1 * (V* 2 -V) t *M y;. 2 -y,')r 
E „ = 5; - P ■ = - (y„'-y,')r 


(D-l) 

(D-2) 

(D-3) 
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= in D . 
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But 


(See Appendix C) 
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Hence: 


M • E 


{(V 2 ~ X i") *i + ( y i*l~ y i") b i} 

{(*.'- v)*. * ( y ;- y /W 

{•i*A«'*i x i'* b j y i*» _b j y /'} 

{ II II * i // 1 // 1 

a.x - a . x. + b. y “b.y. \ 

J n j j j j J \ J 

{ « i * // a n , u it I 

a i x i*2 + b i y i*2-»i y i*i + x )*. y i; 

{ //,!_ II II U ■ II //I 

a.x +b.y -X. y.. +x. + 1 y. > 

J n j J n J^J+l 
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Since 


— li n _ // /> 

Cj - y j + i -x j + i yj 


(See Appendix C) 

= ( 3 j X j" + 2 +b j yj" + 2 +C j) ( a j X n" +b j y n " +C j) < D - U ) 

The above expression can be thought of as the product of two dot 
products: 

b- e * [(-,.*>,• c,) •(*/.,. y/,,. >)][(»,. b >' c i) • «■ <■ *)] < D - 12 > 

This can be written as: 


M ' E n _ ( L j . j + 1 ' P j + 2 ) ( L j . j + 1 ' ^n) 


(D-13) 
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Where 


L . . + j = (a., , c.^ is a three-vector denoting the 

analytical expression of a line which 
is a boundary of the polygonal shadow 
of a surface of the satellite. 

P , „ = /x.'' , y " , 1) is a three-vector denoting the posi- 

j +2 V j +2’ +2’ ) _ „ s r 

tion of the pj +2 vertex. 

Q n s (x n " t y n ", l) is a three -vector denoting the posi- 

tion of the projection of the nth ex- 
periment on the test plane. 


Hence the test can be written as: 


(«» > 0 (D-14) 

for all j=0, 1,2 ... M - 1 . This is identically equal to equation 

(22) in Section IV B. The derivation shown in this appendix proves that 
equation (22) is a valid means of determining whether or not an experi- 
ment is in the shadow of a plane surface in space (i. e. a component 
part of a satellite) providing that the Sun and experiment are on oppo- 
site sides of the plane surface in space. 
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Figure 3— Coordinate System for Solar Array. 



SOLAR ARRAY SHAFT 



Figure 4— Coordinate System for OPEP. 




MAIN BOX SHADOW HISTORY 
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Figure 6— Main Box Shadow History 



EXPERIMENT PACKAGE SHADOW HISTORY 




OPEP SHADOW HISTORY 



Figure 8— OPEP Shadow History 




Figure 10— Solar Heat Inputs to Main Box 



EMITTED HEAT INPUTS TO MAIN BOX 
FOR ONE ORBITAL PERIOD 
ME OF INJECTION INTO ORBIT, 
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SOLAR REFLECTED HEAT INPUTS TO 
BOOM- MOUNTED EXPERIMENT PACKAGES 
FOR ONE ORBITAL PERIOD 
TIME OF INJECTION INTO ORBIT, 

MAY 15, 1964, 8.0 hours U. T. 



EARTH EMITTED HEAT INPUTS TO SOLAR ARRAY 
FOR ONE ORBITAL PERIOD 
TIME OF INJECTION INTO ORBIT: 
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FOR ONE ORBITAL PERIOD 
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EARTH EMITTED HEAT INPUTS TO OPEP 
FOR ONE ORBITAL PERIOD 
TIME OF INJECTION INTO ORBIT, 



Figure 18 — Earth Emitted Heat Inputs to OPEP 




SOLAR REFLECTED HEAT INPUTS TO OPEP 
FOR ONE ORBITAL PERIOD 
TIME OF INJECTION INTO ORBIT, 
MAY 15, 1964, 8.0 hours U. T. 



